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ELECTROMAGNETIC  SCATTERING  PATTERNS  FROM  SINUSOIDAL  SURFACES: 

VERTICAL  POLARIZATION 


INTRODUCTION 

In  this  paper  we  present  an  analysis  and  calculation  of  scattering 
patterns  resulting  from  a  vertically  polarized  beam  of  electromagnetic 
plane  waves  incident  on  a  sinusoidal  surface.  This  model  is  frequently  used 
to  study  scattering  from  periodically  rough  surfaces.  The  horizontally- 
polarized  case  is  treated  in  a  previous  paper  [l}  by  the  authors  which 
contains  more  detailed  discussion  and  references.  The  vector  properties 
of  this  electromagnetic  scattering  problem  are  considered  in  this  communica¬ 
tion  by  extending  the  previous  analysis  to  the  vertically-polarized  case  and 
comparing  the  scattering  patterns  for  the  two  polarizations. 

Electromagnetic  waves  that  are  scattered  from  periodically  rough  surfaces 
have  a  characteristic  dependence  on  the  observation  angle.  In  addition  to 
scattering  at  the  specular  angle,  which  is  characteristic  of  reflection  from 
a  smooth  plane,  these  waves  also  scatter  at  a  discrete  number  of  different 
angles,  due  to  the  rough  nonplanar  nature  of  the  surface. 

The  fields  scattered  from  periodically  rough  surfaces  can  be  computed  by 
employing  the  space-harmonic  representation.  This  representation  expresses 
the  scattered  field  as  a  discrete  sum  of  space  harmonics  (plane  waves) .  The 
scattering  coefficient  associated  with  each  space  harmonic  is  calculated  after 
the  application  of  an  exact  boundary  condition  on  the  surface.  The  polariza¬ 
tion  of  the  incident  electromagnetic  wave  relative  to  the  sinusoidal  surface 

determines  the  type  of  boundary  condition  (Dirichlet  or  Neumann)  that  is  used. 
Note:  Manuscript  submitted  January  15,  1980. 


The  previous  analysis  of  the  scattering  patterns  [1]  has  been  generalized 
to  include  both  vertically  and  horizontally  polarized  electromagnetic  waves. 
The  scattering  of  acoustic  (scalar)  waves  from  a  sinusoidal  surface  has  been 
analyzed  by  DeSanto  [2].  Whitman  and  Schwering  W  have  independently  ob¬ 
tained  a  similar  formulation  for  the  electromagnetic  surface  currents  for 
both  the  horizontal  and  vertical  cases  and  have  compared  their  calculations 
with  the  direct  numerical  solutions  of  the  integral  equations  for  the  sur¬ 
face  currents. 


REPRESENTATION  OF  THE  ELECTROMAGNETIC  FIELDS 


The  electromagnetic  field  scattered  from  a  sinusoidal  surface  is  repre¬ 
sented  by  a  superposition  of  plane  waves  which  consists  of  a  discrete 
spectrum  of  space  harmonics.  In  general,  the  electromagnetic  field  incident 
on  the  surface  can  be  analyzed  in  terms  of  two  orthogonally  polarized  com¬ 
ponents:  horizontal  and  vertical.  In  the  horizontal  case  (denoted  by  a 

(H) 

superscript  (H))  the  total  E  field  is  parallel  to  the  y  axis;  E^  =  ay  '(x,z) 

where  a^  is  a  unit  vector  in  the  y  direction.  In  the  vertical  case  (denoted 

by  a  superscript  (V))  the  total  H  field  is  parallel  to  the  y  axis; 

(V) 

H  «  (x,z).  The  total  field  amplitude  above  the  z  =  0  plane  is  composed 

of  the  incident  and  scattered  field  amplitudes 

ij/q^(x,z)  =  (x,z)  +  (x,z) ,  q£{H,V}.  (1) 


The  spatial  dependence  of  the  incident  plane  wave  is 

ik(a  x-B  z) 


(x , z)  =  e 


z  2  0  (2) 

where  k  ■  ^ /\  is  the  wavenumber  and  a  =  sin  6 .  and  B  **  cos  9 . .  Here  9 . 

o  i  o  i  I 

is  the  angle  of  incidence  measured  with  respect  to  the  z  axis.  The 
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space-harmonic  representation  of  the  scattered  field  is 


*  R  (,)elk(V  +  SnZ),  *  *  0 

sc  n 


where  R  i9  the  scattering  coefficient  of  the  nth  space  harmonic, 
c*n  *  sin0n,  ^n  *  cos®n.>  83111  ®n  1-8  tlie  scatterl-n8  angle  of  the  nth  space  har¬ 
monic  as  shown  in  Figure  1.  We  note  that  the  angle  of  specular  reflection  is 

0  and  0  *  0  . 

o  o  i 

The  profile,  s(x),  of  the  surface  shown  in  Figure  1  is  a  sinusoid  with 
period  l  and  amplitude  d 

s(x)=-^  (1  +  cos  -y-  x).  (4) 

Due  to  the  periodicity  of  the  surface,  the  normalized  transverse  wavenumbers, 
c*n,  satisfy  the  grating  equation. 

o/  *  a  +  nA,  n  »  + 1,  +2,  ....  (5) 

where  A  =  X/i,  and 

6n  =  Jl  "  al  >  Im  (0n}  *  °*  (6) 

The  scattering  coefficient,  Rn^  can  be  calculated  by  applying  Green's  second 

theorem  to  the  region  bounded  by  the  curve  C,  as  is  shown  in  Figure  1.  We 


Jc  [•  •*  i<q)]  *  - ». 


where 


„+  ik(+6  z  -  a  x) 

G—  ■  e  -  n  n 


Here  u  is  the  arc  length  along  the  contour  C  and  n  is  the  inward  normal. 

Along  that  portion  of  the  contour  on  the  sinusoidal  surface,  the  boundary 

condition  is  used  to  simplify  the  integral  in  (7).  The  vanishing  of  the 

tangential  E  field  on  the  surface  of  the  conductor  implies  that  ''j  ’  *  0 

,,  (V) 

(Dirichlet  condition)  and  “■  ■  0  (Neumann  condition)  on  the  surface. 

on 


By  using  the  above  boundary  conditions  in  the  integral  equation  (7)  ex¬ 
pressions  for  the  scattering  coefficients  R  ^  are  obtained  ([3]  and  [2]) 


<1>  ,  ,  .  .  0,  ±1.  ±2, 

m  ^  mn  n  9  9 


(9) 


where 


V  (H)  -  -  im  |m  j  (T  +), 

mn  6  m-n  m 

o 

V  (V)  -  im  (1+b  )  J  (t  +) 

mn  mn  m-n  m 


(10) 

(ID 


a  (a  -  a  ) 
mm  n 

mn  “  6^(6  +  8  ) 

mm  o 


(12) 


Here  is  the  nth-order  Bessel  function  with  argument 


t-  =  A(B  ±  8  ),  A  «  ird/A . 
m  o  m 


(13) 


The  that  appear  in  (9)  are  the  coefficients  of  the  Fourier  expan¬ 

sion  of  the  surface  currents  and  thus  will  be  called  the  surface  current 
amplitudes.  They  satisfy  the  following  infinite  system  of  equations: 

(q)  .  (q) 


T  a  "1'  =  <5  ,  m  *  0,  ±1,  . . . 

mn  n  m,o 


(14) 


where 


T  (H)  =  J  (x“) 
mn  m— n  m 


(15) 


(V) 


mn 


(1  -  8  8 

o  m 


a  a  )J  (t  ) ,  m  =  ±1,  ±2,  ... 
m  n  m-n  m 


(16) 


d  +  c5  ,  +  c<5  .  ,  m  =  0 

n,o  n,l  n,-l 


with  c  =  iAAa  / (28  )  and  6  is  the  Kronecker  delta  function: 
o  o  m,n 


5  =  ii,  m  =  n 

m,n  (0,  m  ¥  n 
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CALCULATION  OF  SCATTERING  COEFFICIENTS  -  VERTICAL  POLARIZATION 


The  infinite  set  of  linear  equations  (14)  can  be  solved  numerically  for 
the  by  symmetrically  truncating  the  system  at  an  order  M  and  solving 

the  resulting  2M-order  matrix  equation.  More  accuracy  is  obtained  by  sequen¬ 
tially  increasing  M;  this  technique  is  known  as  the  method  of  reduction. 

Confidence  in  this  numerical  method  for  solving  the  set  of  equations  (14) 
was  established  by  employing  two  criteria:  conservation  of  energy  and  solu¬ 
tion  stability.  For  a  lossless  surface  the  conservation  of  energy  criterion 

requires  that  the  incident  power  be  equal  to  the  scattered  power.  Employing 

* 

Poynting's  theorem,  we  have  (z  is  the  complex  conjugate  of  z) 


where 

eM  =  I1  "S  irn(M)  !2  I  *  (21) 

n£P 

Here  is  the  energy-conservation  error  made  when  the  approximate  reflection 

coefficient,  T  (M) ,  is  used;  the  r  (M)  are  obtained  when  the  {a  }  have  been 
n  n  n 

computed  after  truncating  system  (14)  to  order  M. 

The  other  criterion  used  was  a  check  of  the  stability  of  the  matrix 
solutions.  The  complex  values  of  a  symmetric  subset  of  the  surface  current 
amplitudes  {a^}  should  stabilize  with  increasing  matrix  size  2M  for  a  given 


set  of  parameters  6^  A,  A.  In  all  cases  which  we  examined,  the  energy- 
conservation  and  current-stabilization  criteria  were  compatible. 

The  above  numerical  procedure  was  carried  out  for  five  different  sets 
of  parameters  in  both  the  horizontally  and  vertically  polarized  cases.  The 
results  are  shown  in  Table  I.  The  first  two  examples  use  the  same  surface 
parameters  that  were  used  in  M.  Example  1.  In  these  cases  both  the  per¬ 
turbation  and  physical-optics  conditions,  as  defined  in  [l] ,  are  exceeded. 

The  last  three  examples  show  Bragg  backscattering.  Examples  three  and  four 
demonstrate  first-order  Bragg  backscattering,  while  example  five  shows 
second-order  Bragg  backscattering.  The  energy-conservation  criterion  was 
applied  with  the  following  steps:  First,  M  was  chosen  large  enough  to 

include  the  set  P  of  all  propagating  modes;  second,  M  was  incremented  until 
-4 

em  <10  ,  if  this  criterion  hadn't  been  satisfied  during  the  calculations 

for  the  first  step.  In  Table  I  the  values  of  normalized  height.  A,  and 
period.  A,  of  the  surface  are  given  along  with  the  angle  of  incidence,  0^. 

We  also  shown  the  maximum  matrix  order  needed  to  establish  an  accuracy  e„. 

Because  of  the  interest  in  radar  backscatter,  we  have  used  our  numerical 
procedure  to  calculate  the  absolute  value  of  the  Bragg  backscattering  ampli¬ 
tude  |r  |.  The  plots  are  shown  in  Figure  2  as  a  function  of  the  normalized 
roughness  amplitude  for  several  different  values  of  Bragg  backscattering 

angle  9  ■  arcsin(A/2).  For  small  values  of  A,  the  backscattering  ampli- 

D,-l 

tude  depends  linearly  on  A:  however,  the  slope  of  this  line  approaches  zero 
as  Qg  -*■  n/2  (grazing).  This  indicates  that  a  perturbation  theory  could 
be  used  for  small  9  .  but  it  would  have  limited  validity  as  9  -*  tt/2.  As 

Dy -1 

A  increases  the  curves  oscillate  between  0  and  1.  The  maximum  value  that  each 
curve  attains  is  1  as  can  be  seen  from  (19).  We  also  note  that  when 
| :_1 I  *  1,  all  other  r  ,  n  ^  -l,  are  zero  because  of  energy  conservation. 
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Thus,  at  this  maximum  all  incident  energy  is  backscattered  to  the  transmitter; 

this  phenomenon  is  known  as  blazing.  It  is  interesting  to  observe  that  all 

curves  shown  in  Figure  2  that  attain  the  value  1  on  the  first  oscillation  have 

only  two  propagating  modes.  Thus  the  curves  for  9_  =  75°  and  85°  do  not 

'  B, -1 

reach  1  on  its  first  oscillation  and  have  more  than  two  modes  propagating. 

We  also  note  that  the  curves  have  been  terminated  at  different  values  of  A; 
this  is  due  to  limitations  of  the  computer  program  whose  convergence  is  para¬ 
meter  dependent. 


CALCULATION  OF  SCATTERING  PATTERNS 


The  scattering  pattern  S^(6)  is  defined  as 

S<q)(0)  =  lim£  |t  (q)|2 

47T  1  SC  1 


(22) 


The  scattered  intensity,  y  4  ,  can  be  calculated  as  a  function  of  observation 

sc 

angle  when  a  bounded  beam  is  incident  on  the  surface.  Let  us  assume  that  the 
incident  beam  is  given  as 


iJ>o(q)  (x,z) 


=  W  jWztane^/Lj 


ik(a  x  -  8  z) 
e  o  o 


(23) 


where  the  beam  is  centered  about  the  line  x  =-z  tan  0^,.  The  beam  shape  is 
defined  by  the  function  W(x/L)  which  is  slowly  varying  when  ]x|<L/2  and 
W(x/L)  =  0  for  |x|>L/2.  If  we  decompose  the  incident  beam  into  plane  waves, 
a  space  harmonic  representation  can  be  obtained  for  the  scattered  field 

(q) 


sc 


An  exact  expression  for  the  scattering  pattern  in  the  horizontally- 
polarized  case  is  derived  In  by  using  the  space  harmonic  representation 
of  the  scattered  field  in  (22).  Since  the  arguments  used  there  hold  equally 
well  in  the  vertically-polarized  case,  we  have 

S^(0)  ■  |^Q  l"kL(sin9-sine  )1  R  ^  (sin8-nA)  |  2 ,  (24) 

n£P  L  n  J  n 


/ 


where 


Q(kLCT)  -  J  W(x)eikLaKdx  .  (25) 

-00 

Here  Q  can  be  interpreted  as  the  scattering  amplitude  due  to  an  incident 
beam  on  the  2=0  plane.  In  addition,  we  note  that  only  propagating  space 
harmonics  contribute  to  S^(0). 

To  use  (24),  one  requires  the  knowledge  of  R  ^  at  many  different 

n 

angles  of  incidence.  When  kL  is  large,  an  approximate  expression  can  be 
obtained  which  only  requires  R  ^  at  the  angle  of  incidence  of  the  beam. 
Following  the  arguments  of  [l],  we  have 

S(q)(0)  =  I Y*  Q(kL(sin0-sin6  ))•  R  (q)(sin6.)|2  (26) 

1  t—J  n  n  i 

n£p 

The  validity  of  (26)  is  based  on  the  fact  that  Rn^(sin0^)  must  be  a  slowly 
varying  function  of  0^  compared  to  the  variation  of  the  Q  function. 

DISCUSSION 

The  exact  formula  (26)  for  the  horizontally-  or  vertically-polarized 

scattering  patterns  will  differ  only  in  the  scattering  coefficients  Rr 

or  R  Thus  the  behavior  of  R  ^  and  R  ^  as  functions  of  the  input 

n  n  n 

parameters  0^,  A,  and  A  (we  recall  that  A  and  A  contain  the  incident  elec¬ 
tromagnetic  wavelength  \)  will  reveal  the  behavior  of  the  corresponding 
scattering  patterns.  For  example,  the  graphs  of  and  as  functions 

of  incidence  angle,  which  we  call  the  space -harmonic  diagrams,  have  been 
extensively  studied  for  periodic  surfaces  (see,  for  example,  [3],  Figs.  3 
and  4).  These  diagrams  reveal  the  existence  and  relative  strengths  of  the 
various  propagating  space  harmonics.  The  graphs  of  R  ^  and  R  ^  as 
functions  of  A  reveal  the  existence  of  blazing;  the  existence  of  Bragg 
backscatter  blazing  is  shown  in  Fig.  2. 
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We  first  repeated  all  the  numerical  examples  of  Table  1  in  [ l]  for  the 
case  of  vertical  polarization.  Since  some  of  these  examples  did  not  differ 
significantly  between  the  two  polarizations,  we  varied  the  input  parameters 
(0^,  A,  A)  to  demonstrate  several  interesting  effects  of  polarization  and 
these  are  summarized  in  Table  1.  However,  in  the  interests  of  brevity,  only 
scattering  patterns  for  the  case  of  blazing  (Example  4  and  Figs.  3  and  4) 
will  be  shown  here. 

Example  1,  provides  a  straightforward  demonstration  of  the  scattering 
pattern  calculation.  However,  for  this  example  the  vertically  and 
horizontally-polarized  cases  happen  to  be  similar  at  9^  =  45°,  which  was 
discussed  in  [l].  Example  1.  An  examination  of  the  space  harmonic  plots 
reveals  that  changing  the  incidence  angle  to  30°,  as  in  example  2,  provides 
a  clearer  distinction  between  the  scattering  patterns  for  the  horizontally- 
and  vertically-polarized  cases. 

The  third  and  fourth  examples  demonstrate  first-order  Bragg  backscatter- 
ing.  The  vertically-polarized  scattering  pattern  for  example  3  is  not 
significantly  different  from  the  corresponding  horizontally-polarized  case 
([l],  Fig.  7)  for  8^  =  45°.  However,  if  the  incidence  angle  is  adjusted  to 
the  blazing  angle,  8^  =  60°,  as  in  example  4,  there  is  substantial  difference 
between  the  horizontally-polarized  (Fig.  3)  and  vertically-polarized  (Fig.  4) 
cases.  Finally,  example  5  demonstrates  the  case  of  vertically-polarized 
second-order  Bragg  backscatter;  comparison  with  the  horizontally-polarized 
case  of  [ l],  Fig.  8,  shows  that  the  specular  and  backscatter  components  are 
larger  for  this  vertical-polarization  case. 
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TABLE  1.  Parameters  for  scattering  pattern  calculations. 


Hor 


Fig.  2.  Bragg  backscatterlng  amplitude,  |r_J  '\,  vs.  roughness  parameter 

for  vertical  polarization. 


8 


6  (DEGREES)-* 

Fig.  4.  Scattering  pattern  for  example  4  (horizontal  polarization). 


